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Abstract 



The origin of cosmic rays is one of the major unresolved astrophysical questions. 
In particular, the highest energy cosmic rays observed possess macroscopic energies 
and their origin is likely to be associated with the most energetic processes in the 
Universe. Their existence triggered a flurry of theoretical explanations ranging from 
conventional shock acceleration to particle physics beyond the Standard Model and 
processes taking place at the earliest moments of our Universe. Furthermore, many new 
experimental activities promise a strong increase of statistics at the highest energies 
and a combination with 7— ray and neutrino astrophysics will put strong constraints 
on these theoretical models. Detailed Monte Carlo simulations indicate that charged 
ultra-high energy cosmic rays can also be used as probes of large scale magnetic fields 
whose origin may open another window into the very early Universe. We give an 
overview over this quickly evolving research field. 
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1 Introduction 
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Fi gure 1: The cosmic ray all particle spectrum Q. Approximate integral fluxes are also shown 



After almost 90 years of research on cosmic rays (CRs), their origin is still an open 
question, for which the degree of uncertainty increases with energy: Only below 100 MeV 
kinetic energy, where the solar wind shields protons coming from outside the solar system, 
the sun must give rise to the observed proton flux. The bulk of the CRs up to at least an 
energy of = 4 x 10^^ eV is believed to originate within our Galaxy. Above that energy, 
which is associated with the so called "knee", the flux of particles per area, time, solid 
angle, and energy, which can be well approximated by broken power laws oc , steepens 
from a power law index 7 ^ 2.7 to one of index ~ 3.2. Above the so called "ankle" at 
~ 5 X 10^^ eV, the spectrum flattens again to a power law of index 7 ~ 2.8. This latter 
feature is often interpreted as a cross over from a steeper Galactic component to a harder 
component of extragalactic origin. Fig. |l] shows the measured CR spectrum above 100 MeV, 
up to 3 X 10^° eV, the highest energy measured so far for an individual CR. 

The conventional scenario assumes that all high energy charged particles are accelerated 
in magnetized astrophysical shocks, whose size and typical magnetic field strength determines 
the maximal achievable energy, similar to the situation in man made particle accelerators. 
The most likely astrophysical accelerators for CR up to the knee, and possibly up to the ankle 
are the shocks associated with remnants of past Galactic supernova explosions, whereas for 
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the presumed extragalactic component powerful objects such as active galactic nuclei are 
envisaged. 
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Figure 2: The cosmic ray spectrum above 10^''' eV (from Ref [|T^).The "ankle" is visible at i? ~ 5 x 10^^ eV. 

The main focus of this contribution will be on ultrahigh energy cosmic rays (UHECRs), 
those with energy > 10^^ eV 0, ^, ^ |^ ||, see Fig. ^ For more details on CRs in general 
the reader is referred to recent monographs |Tl| , |12|. In particular, extremely high energy 
(EHE)[] cosmic rays pose a serious challenge for conventional theories of CR origin based on 
acceleration of charged particles in powerful astrophysical objects. The question of the origin 



^We shall use the abbreviation EHE to specifically denote energies E > f 0^° eV, while the abbreviation 
UHE for "Ultra-High Energy" will sometimes be used to denote E > 1 EeV, where 1 EeV = 10^* eV. Clearly 
UHE includes EHE but not vice versa. 



3 



of these EHECRs is, therefore, currently a subject of much intense debate and discussions 
as well as experimental efforts; see Refs. P, ^ |13|, and Ref. |10|, |14| for recent brief reviews, 
and Ref. for a detailed review. In Sect. 2 we will summarize detection techniques and 
present and future experimental projects. 

The current theories of origin of EHECRs can be broadly categorized into two distinct 
"scenarios": the "bottom-up" acceleration scenario, and the "top-down" decay scenario, 
with various different models within each scenario. As the names suggest, the two scenarios 
are in a sense exact opposite of each other. The bottom-up scenario is just an extension of 
the conventional shock acceleration scenario in which charged particles are accelerated from 
lower energies to the requisite high energies in certain special astrophysical environments. On 
the other hand, in the top-down scenario, the energetic particles arise simply from decay of 
certain sufficiently massive particles originating from physical processes in the early Universe, 
and no acceleration mechanism is needed. 

The problems encountered in trying to explain EHECRs in terms of acceleration mech- 



anisms have been well-documented in a number of studies; see, e.g., Refs. [0, |T^, |T8[. Even 
if it is possible, in principle, to accelerate particles to EHECR energies of order 100 EeV 
in some astrophysical sources, it is generally extremely difficult in most cases to get the 
particles come out of the dense regions in and/or around the sources without losing much 
energy. Currently, the most favorable sources in this regard are perhaps a class of powerful 
radio galaxies (see, e.g., Refs. |T^, |2D[ for recent reviews and references to the literature), 
although the values of the relevant parameters required for acceleration to energies > 100 
EeV are somewhat on the extreme side [1^. However, even if the requirements of energetics 



are met, the main problem with radio galaxies as sources of EHECRs is that most of them 
seem to lie at large cosmological distances, ^ 100 Mpc, from Earth. This is a major prob- 
lem if EHECR particles are conventional particles such as nucleons or heavy nuclei. The 
reason is that nucleons above ^ 70 EeV lose energy drastically during their propagation 



from the source to Earth due to the Greisen-Zatsepin-Kuzmin (GZK) effect [21, 22], namely, 



photo-production of pions when the nucleons collide with photons of the cosmic microwave 
background (CMB), the mean- free path for which is ~ few Mpc |23|. This process limits the 
possible distance of any source of EHE nucleons to < 100 Mpc. If the particles were heavy 
nuclei, they would be photo-disintegrated p^ , p5[] in the CMB and infrared (IR) background 
within similar distances. Thus, nucleons or heavy nuclei originating in distant radio galaxies 
are unlikely to survive with EHECR energies at Earth with any significant flux, even if they 
were accelerated to energies of order 100 EeV at source. In addition, if cosmic magnetic 
fields are not close to existing upper limits, EHECRs are not likely to be deflected strongly 
by large scale cosmological and/or Galactic magnetic fields Thus, EHECR arrival directions 
should point back to their sources in the sky (see Sect. 5 for details) and EHECRs may 
offer us the unique opportunity of doing charged particle astronomy. Yet, for the observed 
EHECR events so far, no powerful sources close to the arrival directions of individual events 
are found within about 100 Mpc [T^. Very recently, it has been suggested by Boldt and 



Ghosh that particles may be accelerated to energies ~ 10^^ eV near the event horizons of 
spinning supermassive black holes associated with presently inactive quasar remnants whose 
numbers within the local cosmological Universe (i.e., within a GZK distance of order 50 
Mpc) may be sufficient to explain the observed EHECR flux. This would solve the problem 
of absence of suitable currently active sources associated with EHECRs. A detailed model 
incorporating this suggestion, however, remains to be worked out. 
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There are, of course, ways to avoid the distance restriction imposed by the GZK effect, 
provided the problem of energetics is somehow solved separately and provided one allows new 
physics beyond the Standard Model of particle physics; we shall discuss those suggestions in 
Sect. 3. 

On the other hand, in the top-down scenario, which will be discussed in Sect. 4, the 
problem of energetics is trivially solved from the beginning. Here, the EHECR particles owe 
their origin to decay of some supermassive "X" particles of mass nix ^ 10^° eV, so that 
their decay products, envisaged as the EHECR particles, can have energies all the way up to 
~ nix- Thus, no acceleration mechanism is needed. The sources of the massive X particles 
could be topological defects such as cosmic strings or magnetic monopoles that could be 
produced in the early Universe during symmetry-breaking phase transitions envisaged in 
Grand Unified Theories (GUTs). In an inflationary early Universe, the relevant topological 
defects could be formed at a phase transition at the end of inflation. Alternatively, the X 
particles could be certain supermassive metastable relic particles of lifetime comparable to or 
larger than the age of the Universe, which could be produced in the early Universe through, 
for example, particle production processes associated with inflation. Absence of nearby 
powerful astrophysical objects such as AGNs or radio galaxies is not a problem in the top- 
down scenario because the X particles or their sources need not necessarily be associated 
with any speciflc active astrophysical objects. In certain models, the X particles themselves 
or their sources may be clustered in galactic halos, in which case the dominant contribution 
to the EHECRs observed at Earth would come from the X particles clustered within our 
Galactic Halo, for which the GZK restriction on source distance would be of no concern. 

By focusing primarily on "non-conventional" scenarios involving new particle physics be- 
yond the electroweak scale, we do not wish to give the wrong impression that these scenarios 
explain all aspects of EHECRs. In fact, as we shall see below, essentially each of the spe- 
ciflc models that have been studied so far has its own peculiar set of problems. Indeed, the 
main problem of non-astrophysical solutions of the EHECR problem in general is that they 
are highly model dependent. On the other hand, it is precisely because of this reason that 
these scenarios are also attractive — they bring in ideas of new physics beyond the Standard 
Model of particle physics (such as Grand Uniflcation and new interactions beyond the reach 
of terrestrial accelerators) as well as ideas of early Universe cosmology (such as topological 
defects and/or massive particle production in inflation) into the realms of EHECRs where 
these ideas have the potential to be tested by future EHECR experiments. 

The physics and astrophysics of UHECRs are intimately linked with the emerging fleld 
of neutrino astronomy (for reviews see Refs. p8| , p9[] ) as well as with the already established 



fleld of 7— ray astronomy (for reviews see, e.g., Ref. ||3^) which in turn are important sub- 
disciplines of particle astrophysics (for a review see, e.g., Ref. f^). Indeed, as we shall 



see, all scenarios of UHECR origin, including the top-down models, are severely constrained 
by neutrino and 7— ray observations and limits. In turn, this linkage has important con- 
sequences for theoretical predictions of fluxes of extragalactic neutrinos above a TeV or so 
whose detection is a major goal of next-generation neutrino telescopes (see Sect. 2): If these 
neutrinos are produced as secondaries of protons accelerated in astrophysical sources and 
if these protons are not absorbed in the sources, but rather contribute to the UHECR flux 
observed, then the energy content in the neutrino flux can not be higher than the one in 



UHECRs, leading to the so called Waxman Bahcall bound ||3^, |3^ . If one of these assump 



tions does not apply, such as for acceleration sources that are opaque to nucleons or in the 
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TD scenarios where X particle decays produce much fewer nucleons than 7— rays and neu- 
trinos, the Waxman Bahcall bound does not apply, but the neutrino flux is still constrained 
by the observed diffuse 7— ray flux in the GeV range (see Sect. 4.4). 

Finally, in Sect. 5 we shall discuss how, apart from the unsolved problem of the source 
mechanism, EHECR observations have the potential to yield important information on 
Galactic and extragalactic magnetic flelds. 



2 Present and Future UHE CR and Neutrino Experi- 
ments 

The CR primaries are shielded by the Earth's atmosphere and near the ground reveal their 
existence only by indirect effects such as ionization. Indeed, it was the height dependence 
of this latter effect which lead to the discovery of CRs by Hess in 1912. Direct observation 
of CR primaries is only possible from space by flying detectors with balloons or spacecraft. 
Naturally, such detectors are very limited in size and because the differential CR spectrum 
is a steeply falling function of energy (see Fig. |I]), direct observations run out of statistics 
typically around a few 100 TeV. 

Above ~ 100 TeV, the showers of secondary particles created in the interactions of the 
primary CR with the atmosphere are extensive enough to be detectable from the ground. In 
the most traditional technique, charged hadronic particles, as well as electrons and muons 



in these Extensive Air Showers (EAS) are recorded on the ground with standard instru- 
ments such as water Cherenkov detectors used in the old Volcano Ranch |^ and Haverah 
Park experiments, and scintillation detectors which are used now-a-days. Currently oper- 
ating ground arrays for UHECR EAS are the Yakutsk experiment in Russia and the Akeno 
Giant Air Shower Array (AGASA) near Tokyo, Japan, which is the largest one, covering an 
area of roughly 100 km^ with about 100 detectors mutually separated by about 1km 0. 
The Sydney University Giant Air Shower Recorder (SUGAR) ||^ operated until 1979 and 
was the largest array in the Southern hemisphere. The ground array technique allows one 
to measure a lateral cross section of the shower proflle. The energy of the shower-initiating 
primary particle is estimated by appropriately parametrizing it in terms of a measurable 
parameter; traditionally this parameter is taken to be the particle density at 600 m from 
the shower core, which is found to be quite insensitive to the primary composition and the 
interaction model used to simulate air showers. 

The detection of secondary photons from EAS represents a complementary technique. 
The experimentally most important light sources are the fluorescence of air nitrogen excited 
by the charged particles in the EAS and the Cherenkov radiation from the charged particles 
that travel faster than the speed of light in the atmospheric medium. The flrst source is 
practically isotropic whereas the second one produces light strongly concentrated on the 
surface of a cone around the propagation direction of the charged source. The fluorescence 
technique can be used equally well for both charged and neutral primaries and was flrst used 
by the Fly's Eye detector |§] and will be part of several future projects on UHECRs (see 
below). The primary energy can be estimated from the total fluorescence yield. Information 
on the primary composition is contained in the column depth X^ax (measured in gcm~^) at 
which the shower reaches maximal particle density. The average of Xmax is related to the 
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primary energy E by 



Here, Xq is called the elongation rate and Eq is a characteristic energy that depends on 
the primary composition. Therefore, if X^ax and Xq are determined from the longitudinal 
shower profile measured by the fluorescence detector, then Eq and thus the composition, can 
be extracted after determining the energy E from the total fluorescence yield. Comparison 
of CR spectra measured with the ground array and the fluorescence technique indicate 
systematic errors in energy calibration that are generally smaller than ~ 40%. For a more 
detailed discussion of experimental EAS analysis with the ground array and the fluorescence 
technique see, e.g., Refs. p6| . 



As an upscaled version of the old Fly's Eye Cosmic Ray experiment, the High Resolution 
Fly's Eye detector is currently under construction at Utah, USA [0. Taking into account 
a duty cycle of about 10% (a fluorescence detector requires clear, moonless nights), the 
effective aperture of this instrument will be ~ 350(1000) km^ sr at 10(100) EeV, on average 
about 6 times the Fly's Eye aperture, with a threshold around 10^^ eV. Another project 
utilizing the fluorescence technique is the Japanese Telescope Array |^ which is currently 
in the proposal stage. If approved, its effective aperture will be about 10 times that of 
Fly's Eye above 10^^ eV, and it would also be used as a Cherenkov detector for TeV 7— ray 
astrophysics. The largest project presently under construction is the Pierre Auger Giant 



Array Observatory planned for two sites, one in Argentina and another in the USA 
for maximal sky coverage. Each site will have a 3000 km^ ground array. The southern 
site will have about 1600 particle detectors (separated by 1.5 km each) overlooked by four 
fluorescence detectors. The ground arrays will have a duty cycle of nearly 100%, leading 
to an effective aperture about 30 times as large as the AGASA array. The corresponding 
cosmic ray event rate above 10^° eV will be about 50 events per year. About 10% of the 
events will be detected by both the ground array and the fluorescence component and can 
be used for cross calibration and detailed EAS studies. The energy threshold will be around 
10^^ eV, with full sensitivity above 10^^ eV. 

Recently NASA initiated a concept study for detecting EAS from space [^, ^ by ob- 
serving their fluorescence light from an Orbiting Wide-angle Light-collector (OWL). This 
would provide an increase by another factor ~ 50 in aperture compared to the Pierre Auger 
Project, corresponding to an event rate of up to a few thousand events per year above 10^° eV. 
Similar concepts such as the Extreme Universe Space Observatory (EUSO) which is part 



of the AirWatch program |^5[ and of which a prototype may be tested on the International 
Space Station are also being discussed. It is possible that the OWL and AirWatch efforts will 
merge. The energy threshold of such instruments would be between 10^^ and 10^" eV. This 
technique would be especially suitable for detection of very small event rates such as those 
caused by UHE neutrinos which would produce deeply penetrating EAS (see Sect. 4.4). For 
more details on these recent experimental considerations see Ref. ITB . 



High energy neutrino astronomy is aiming towards a kilometer scale neutrino observatory. 
The major technique is the optical detection of Cherenkov light emitted by muons created 
in charged current reactions of neutrinos with nucleons either in water or in ice. The largest 
pilot experiments representing these two detector media are the now defunct Deep Undersea 
Muon and Neutrino Detection (DUMAND) experiment |4^ in the deep sea near Hawai and 
the Antarctic Muon And Neutrino Detector Array (AMANDA) experiment in the South 
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Pole ice. Another water based experiment is situated at Lake Baikal |4^. Next generation 
deep sea projects include the French Astronomy with a Neutrino Telescope and Abyss envi- 
ronmental RESearch (ANTARES) ||4^ and the underwater Neutrino Experiment SouthwesT 
Of GReece (NESTOR) project in the Mediterranean [Q, whereas ICECUBE represents 
the planned kilometer scale version of the AMANDA detector. Also under consideration are 
neutrino detectors utilizing techniques to detect the radio pulse from the electromagnetic 
showers created by neutrino interactions in ice. This technique could possibly be scaled up 
to an effective area of 10^ km^ and a prototype is represented by the Radio Ice Cherenkov 
Experiment (RICE) experiment at the South Pole |52|. Neutrinos can also initiate horizontal 
EAS which can be detected by giant ground arrays such as the Pierre Auger Project [p3|. 



Furthermore, as mentioned above, deeply penetrating EAS could be detected from space by 
instruments such as the proposed space based Air Watch type detectors . More 

details and references on neutrino astronomy detectors are contained in Refs. PSI, and 



some recent overviews on neutrino astronomy can be found in Ref. p9 



3 New Primary Particles and New Interactions 

A possible way around the problem of missing counterparts within acceleration scenarios 
is to propose primary particles whose range is not limited by interactions with the CMB. 
Within the Standard Model the only candidate is the neutrino, whereas in supersymmetric 
extensions of the Standard Model, new neutral hadronic bound states of light gluinos with 
quarks and gluons, so-called R-hadrons that are heavier than nucleons, and therefore have 



a higher GZK threshold, have been suggested |55 



In both the neutrino and new massive neutral hadron scenario the particle propagating 
over extragalactic distances would have to be produced as a secondary in interactions of a 
primary proton that is accelerated in a powerful AGN which can, in contrast to the case 
of EAS induced by nucleons, nuclei, or 7— rays, be located at high redshift. Consequently, 
these scenarios predict a correlation between primary arrival directions and high redshift 
sources. In fact, possible evidence for an angular correlation of the five highest energy events 
with compact radio quasars at redshifts between 0.3 and 2.2 was recently reported [^]. A 



new analysis with the somewhat larger data set now available does not support significant 
correlations This is currently disputed since another group claims to have found a 

correlation on the 99.9% confidence level [Q. Only a few more events could confirm or rule 
out the correlation hypothesis. Note that these scenarios would require the primary proton 
to be accelerated up to > 10^^ eV, demanding a very powerful astrophysical accelerator. On 
the other hand, a few dozen such exceptional accelerators in the visible Universe may suffice. 



3.1 New Neutrino Interactions 

Neutrino primaries have the advantage of being well established particles. However, within 
the Standard Model their interaction cross section with nucleons, whose charged current part 
can be parametrized by |^ 

a^J{E) ^ 2.36 X IQ-^^E/IQ^'^ ^Yf-^^^ cm^ (10^^ eV < E < 10^^ eV) , (2) 
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falls short by about five orders of magnitude to produce ordinary air showers. However, it 
has been suggested that the neutrino-nucleon cross section, (Tj^at, can be enhanced by new 
physics beyond the electroweak scale in the center of mass (CM) frame, or above about a 
PeV in the nucleon rest frame. Neutrino induced air showers may therefore rather directly 
probe new physics beyond the electroweak scale. 

Two major possibilities have been discussed in the literature for which unitarity bounds 
need not be violated. In the first, a broken SU(3) gauge symmetry dual to the unbroken 
SU(3) color gauge group of strong interaction is introduced as the "generation symmetry" 
such that the three generations of leptons and quarks represent the quantum numbers of this 
generation symmetry. In this scheme, neutrinos can have close to strong interaction cross 
sections with quarks. In addition, neutrinos can interact coherently with all partons in the 
nucleon, resulting in an effective cross section comparable to the geometrical nucleon cross 
section. This model lends itself to experimental verification through shower development 
altitude statistics [|60| . 

The second possibility consists of a large increase in the number of degrees of freedom 
above the electroweak scale pT|. A specific implementation of this idea is given in theories 



with n additional large compact dimensions and a quantum gravity scale M^+n ~ TeV that 



has recently received much attention in the literature ||62[ because it provides an alternative 
solution (i.e., without supersymmetry) to the hierarchy problem in grand unifications of 
gauge interactions. The cross sections within such scenarios have not been calculated from 
first principles yet. Within the field theory approximation which should hold for squared 
CM energies s < M|_,_„, the spin 2 character of the graviton predicts ag ~ s'^/Mf_^^ |^ For 
s ^ ^4+„, several arguments based on unitarity within field theory have been put forward. 



Ref. [53 1 suggested 



where in the last expression we specified to a neutrino of energy E hitting a nucleon at rest. 
A more detailed calculation taking into account scattering on individual partons leads to 



similar orders of magnitude ||6^. Note that a neutrino would typically start to interact in 
the atmosphere for a^N ^ 10~^^cm^, i.e. in the case of Eq. for E > 10^° eV, assum- 
ing M^^n — 1 TeV. The neutrino therefore becomes a primary candidate for the observed 
EHECR events. However, since in a neutral current interaction the neutrino transfers only 
about 10% of its energy to the shower, the cross section probably has to be at least a few 
10^^^ cm^ to be consistent with observed showers which start within the first 50gcm~^ of 



the atmosphere p3[. A specific signature of this scenario would be the absence of any events 



above the energy where ag grows beyond ~ 10 ^^cm^ in neutrino telescopes based on ice 
or water as detector medium [2I5|, and a hardening of the spectrum above this energy in 



atmospheric detectors such as the Pierre Auger Project |^ and the proposed space based 
Air Watch type detectors |^2|, ^ 0, Furthermore, according to Eq. (^, the average 
atmospheric column depth of the first interaction point of neutrino induced EAS in this 
scenario is predicted to depend linearly on energy. This should be easy to distinguish from 
the logarithmic scaling, Eq. (|I]), expected for nucleons, nuclei, and 7— rays. To test such 
scalings one can, for example, take advantage of the fact that the atmosphere provides a 
detector medium whose column depth increases from ~ 1000 g/cm^ towards the zenith to 
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36000 g/cm^ towards horizontal arrival directions. This probes cross sections in the range 
X0~29 _ io~27j.j^2_ -Q^g ^Yie increased Water/ice detectors would probe cross sections 

6§. 



-31 



in the range ~ 10 



10 



-29 



Within string theory, individual amplitudes are expected to be suppressed exponentially 
above the string scale Ms which for simplicity we assume here to be comparable to M4+„. 
This can be interpreted as a result of the finite spatial extension of the string states. In 
this case, the neutrino nucleon cross section would be dominated by interactions with the 
partons carrying a momentum fraction x ~ M^/s, leading to 



a.N ^ —Hs/M^){s/Ml 



6 X 10-2^ 



2^0.363 



-4.726 



E 



X 



1 + 0.08 In 



TeVy VlO^OeV 
E 



0.363 



(4) 



1020 eV 



0.16 In 



Ms \ 
TeVy 



cm 



This is probably too small to make neutrinos primary candidates for the highest energy 
showers observed, given the fact that complementary constraints from accelerator experi- 
ments result in Mg > 1, TeV |^|. On the other hand, in the total cross section amplitude 
suppression may be compensated by an exponential growth of the level density |^1|. It is 
currently unclear and it may be model dependent which effect dominates. Thus, an experi- 
mental detection of the signatures discussed in this section could lead to constraints on some 
string-inspired models of extra dimensions. 

We note in passing that extra dimensions can have other astrophysical ramifications such 
as energy loss in stellar environments due to emission of real gravitons into the bulk. The 
strongest resulting lower limits on M4+„ come from the consideration of cooling of the cores 
of hot supernovae and read Mg > 50 TeV, M7 > 4 TeV, Mg > 1 TeV, and Mu > 0.05 TeV 
for n = 2,3,4,7, respectively In addition, implications of extra dimensions for early 

Universe physics and inflation are increasingly studied in the literature, but much work is 
left to be done on the intersection of these research domains. 

Independent of theoretical arguments, the UHECR data can be used to put constraints 
on cross sections satisfying a^N^E > 10^^ eV) < lO'^'^cm^. Particles with such cross sections 
would give rise to horizontal air showers. The Fly's Eye experiment established an upper 
limit on horizontal air showers . The non-observation of the neutrino flux expected from 



pions produced by UHECR interacting with the CMB the results in the limit 



aMlO^'eV) < 8x10-^7^^2 
a.^(lO^^V) < 5 X 10-^7y ^2 cm^ 



cm 



cm 



(5) 



where y is the average energy fraction of the neutrino deposited into the shower {y = 1 for 
charged current reactions and ^ ~ 0.1 for neutral current reactions). Expected neutrino 
fluxes are shown in Fig. ^. The projected sensitivity of future experiments such as the Pierre 
Auger Observatories and the AirWatch type satellite projects indicate that the cross section 
limits Eq. (|) could be improved by up to four orders of magnitude, corresponding to one 
order of magnitude in Mg or M4+„. 
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3.2 Supersymmetric Particles 

Light gluinos binding to quarks, anti-quarks and/or gluons can occur in supersymmetric 



theories involving gauge-mediated supersymmetry (SUSY) breaking [|7T| where the resulting 
gluino mass arises dominantly from radiative corrections and can vary between ~ 1 GeV 
and ~ 100 GeV. In these scenarios, the gluino can be the lightest supersymmetric particle 
(LSP). There are also arguments against a light quasi-stable gluino [^, mainly based on 
constraints on the abundance of anomalous heavy isotopes of hydrogen and oxygen which 
could be formed as bound states of these nuclei and the gluino. Furthermore, accelerator 
constraints have become quite stringent 17^ and seem to be inconsistent with the original 



scenario from Ref. ||5^. However, the scenario with a "tunable" gluino mass |^ still seems 
possible and suggests either the gluino-gluon bound state gg, called glueballino Rq, or the 
isotriplet g — {uu — dd)s, called p, as the lightest quasi-stable R-hadron. For a summary of 



scenarios with light gluinos consistent with accelerator constraints see Ref. The case of 
a light quasi-stable gluino does not seem to be settled. 

An astrophysical constraint on new neutral massive and strongly interacting EAS pri- 
maries results from the fact that the nucleon interactions producing these particles in the 
source also produce neutrinos and especially 7— rays. The resulting fluxes from powerful dis- 
crete acceleration sources may be easily detectable in the GeV range by space-borne 7— ray 
instruments such as EGRET and GLAST, and in the TeV range by ground based 7— ray 
detectors such as HEGRA and WHIPPLE and the planned VERITAS, HESS, and MAGIC 
projects (for reviews discussing these instruments see Ref. fSyj). At least the latter three 
ground based instruments should have energy thresholds low enough to detect 7— rays from 
the postulated sources at redshift z ~ 1. Such observations in turn imply constraints on 
the required branching ratio of proton interactions into the R-hadron which, very roughly, 
should be larger than ~ 0.01. These constraints, however, will have to be investigated in 
more detail for specific sources. One could also search for heavy neutral baryons in the 
data from Cherenkov instruments in the TeV range in this context. To demonstrate these 
points, a schematic example of fluxes predicted for the new heavy particle and for 7— rays 
and neutrinos are shown in Fig. ^ 

A further constraint on new EAS primary particles in general comes from the character 
of the air showers created by them: The observed EHECR air showers are consistent with 
nucleon primaries and limits the possible primary rest mass to less than ~ 50 GeV . With 
the statistics expected from upcoming experiments such as the Pierre Auger Project, this 
upper limit is likely to be lowered down to ~ 10 GeV. 

It is interesting to note in this context that in case of a confirmation of the existence of 
new neutral particles in UHECRs, a combination of accelerator, air shower, and astrophysics 
data would be highly restrictive in terms of the underlying physics: In the above scenario, 
for example, the gluino would have to be in a narrow mass range, 1-10 GeV, and the newest 
accelerator constraints on the Higgs mass, ruh > 90 GeV, would require the presence of a D 
term of an anomalous U{l)x gauge symmetry, in addition to a gauge-mediated contribution 



to SUSY breaking at the messenger scale . 
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Figure 3: Schematic predictions for the fluxes of the putative new neutral heavy particle (dotted Une), 
electron, muon, and r— neutrinos (dashed and dash-dotted lines, as indicated), and 7— rays (solid line) for 
a source at redshift z = 1. Assumed were a proton spectrum cx E^"^-^ extending at least up to 10^^ eV 
at the source, a branching ratio for production of the heavy neutral in nucleon interactions of 0.01, and a 
beaming factor of 10 for neutrinos and the heavy neutrals. The 1 sigma error bar at 3 x 10^" eV represents 
the point flux corresponding to the highest energy Fly's Eye event. The predicted fluxes were normalized 
such that this highest energy event is explained as a new heavy particle. The points with arrows on the right 
part represent projected approximate neutrino point source sensitivities for the space based AirWatch type 
concepts using the OWL/ Air Watch acceptance estimated in Ref. |^2[ ^ for non-detection over a five year 
period. The points with arrows in the lower left part represent approximate 7— ray point source sensitivities 
of existing detectors such as EGRET and HEGRA, and of planned instruments such as the satellite detector 
GLAST, the Cherenkov telescope array HESS and the single dish instrument MAGIC, for 50 hours and 1 
month observation time for the ground based and satellite detectors, respectively. 



3.3 Anomalous Kinematics, Quantum Gravity Effects, Lorentz 
Symmetry Violations 

The existence of UHECR beyond the GZK cutoff has prompted several suggestions of pos- 
sible new physics beyond the Standard Model. We have already discussed some of these 
suggestions in Sect. 3.1 and 3.2 in the context of new primary particles and new interactions. 
Further, in Sect. 4 we will discuss suggestions regarding possible new sources of EHECR that 
also involve postulating new physics beyond the Standard Model. In the present section, to 
end our discussions on the propagation and new interactions of UHE radiation, we briefly 
discuss some examples of possible small violations or modifications of certain fundamental 
tenets of physics (and constraints on the magnitude of those violations/modifications) that 
have also been discussed in the literature in the context of propagation of UHECR. 

For example, as an interesting consequence of the very existence of UHECR, constraints 



on possible violations of Lorentz invariance (VLI) have been pointed out [Q. These con 
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straints rival precision measurements in the laboratory: If events observed around 10^° eV 
are indeed protons, then the difference between the maximum attainable proton velocity 
and the speed of light has to be less than about 1 x 10"^'^, otherwise the proton would lose 
its energy by Cherenkov radiation within a few hundred centimeters. Possible tests of other 
modes of VLI with UHECR have been discussed in Ref. IfrTl, and in Ref. HTSI in the context of 



horizontal air-showers generated by cosmic rays in general. Gonzalez-Mestres [77|, Coleman 



and Glashow [^], and earlier, Sato and Tati |^ and Kirzhnits and Chechin have also 



suggested that due to modified kinematical constraints the GZK cutoff could even be evaded 
by allowing a tiny VLI too small to have been detected otherwise. Similar consequences 
apply to other energy loss processes such as pair production by photons above a TeV with 
the low energy photon background [^. It seems to be possible to accomodate such effects 
within theories involving generalized Lorentz transformations |^| or deformed relativistic 
kinematics IS^. Furthermore, it has been pointed out that violations of the principle 



of equivalence (VPE), while not dynamically equivalent, also produce the same kinematical 
effects as VLI for particle processes in a constant gravitational potential, and so the con- 
straints on VLI from UHECR physics can be translated into constraints on VPE such that 
the difference between the couplings of protons and photons to gravity must be less than 
about 1 X 10^^^. Again, this constraint is more stringent by several orders of magnitude 
than the currently available laboratory constraint from Eotvos experiments. 

As a specific example of VLI, we consider an energy dependent photon group velocity 
dE/dk = c[l — xE/Eq + OlE"^ / Eq)] where c is the speed of light in the low energy limit, 
X = ±1, and Eq denotes the energy scale where this modification becomes of order unity. 
This corresponds to a dispersion relation 



2u2 



^2 



(6) 



which, for example, can occur in quantum gravity and string theory |86|. The kinematics of 
electron-positron pair production in a head-on collision of a high energy photon of energy E 
with a low energy background photon of energy e then leads to the constraint 



EiE"^ 



+ 9102 



X 



E^ 

4:Ea 



(7) 



where Ei and 6i ~ 0{m/Ei) are respectively the energy and outgoing momentum angle 
(with respect to the original photon momentum) of the electron and positron {i = 1, 2). For 
the case considered by Coleman and Glashow |7B| in which the maximum attainable speed 
Ci of the matter particle is different from the photon speed c, the kinematics can be obtained 
by substituting cf — for xE/Eq in Eq. (|^. 

Let us define a critical energy E^ = {mlEQy^^ ^ Ib^Eo/mpi)^/^ TeV in the case of the 
energy dependent photon group velocity, and Ec = rne/\c^ — c^|^/^ in the case considered 
by Coleman and Glashow. If x < 0, or Cj < c, then e becomes negative for E > E^.. 
This signals that the photon can spontaneously decay into an electron-positron pair and 
propagation of photons across extragalactic distances will in general be inhibited. The 
observation of extragalactic photons up to ~ 20 TeV ^] therefore puts the limits Eq > 



Mpi or cf — > — 2 X 10 In contrast, if x > 0, or Cj > c, e will grow with energy 
for E > Ec until there is no significant number of target photon density available and the 



13 



Universe becomes transparent to UHE photons. A clear test of this possibihty would be the 
observation of > 100 TeV photons from distances > 100 Mps [^. An unambiguous detection 
of the thresholds for pion production by nucleons and for pair production by photons in the 
CMB and other low energy photon backgrounds in the future would allow to establish more 
stringent lower limits on Eq |]90| . 

In addition, the dispersion relation Eq. (^) implies that a photon signal at energy E will 
be spread out by At ~ {d/c){E/Eo) ^ l(rf/100 Mpc)(E/ TeV)(Eo/Mpi)-i s. The observation 
of 7— rays at energies E > 2 TeV within ~ 300 s from the AGN Markarian 421 therefore puts 
a limit (independent of x) of Eq > Ax 10^^ GeV, whereas the possible observation of 7— rays 
atE> 200 TeV within ~ 200 s from a ORB by HEGRA might be sensitive to Eq ^ Mpi [|T |. 



For a recent detailed discussion of these limits see Ref. ^2 . 

A related proposal originally due to Kostelecky in the context of CR suggests the electron 
neutrino to be a tachyon This would allow the proton in a nucleus of mass m{A, Z) 
for mass number A and charge Z to decay via p — n + e"*" + z/g above the energy threshold 
Eth = m{A, Z)[m{A, Z ±1) + mg — m{A, Z)]/\m^J which, for a free proton, is E^h — 1-7 x 
10^^/(|m,^^|/eV) eV. Ehrlich claims that by choosing m^^ ~ — (0.5eV)^ it is possible to 
explain the knee and several other features of the observed CR spectrum, including the high 
energy end, if certain assumptions about the source distribution are made. The experimental 
best fit values of m^^ from tritium beta decay experiments are indeed negative , although 
this is most likely due to unresolved experimental issues. In addition, the values of \ml^ \ from 
tritium beta decay experiments are typically larger than the value required to fit the knee 



of the CR spectrum. This scenario also predicts a neutron line around the knee energy [96 



4 Top-Down Scenarios 
4.1 The Main Idea 

As mentioned in the introduction, all top-down scenarios involve the decay of X particles of 
mass close to the GUT scale which can basically be produced in two ways: If they are very 
short hved, as usually expected in many GUTs, they have to be produced continuously. The 
only way this can be achieved is by emission from topological defects left over from cosmo- 
logical phase transitions that may have occurred in the early Universe at temperatures close 
to the GUT scale, possibly during reheating after inflation. Topological defects necessarily 
occur between regions that are causally disconnected, such that the orientation of the or- 
der parameter associated with the phase transition can not be communicated between these 
regions and consequently will adopt different values. Examples are cosmic strings (similar 
to vortices in superfluid helium), magnetic monopoles, and domain walls (similar to Bloch 
walls separating regions of different magnetization in a ferromagnet). The defect density is 
consequently given by the particle horizon in the early Universe and their formation can even 
be studied in solid state experiments where the expansion rate of the Universe corresponds 



to the quenching speed with which the phase transition is induced |^^. The defects are 
topologically stable, but in the cosmological case time dependent motion leads to the emis- 
sion of particles with a mass comparable to the temperature at which the phase transition 
took place. The associated phase transition can also occur during reheating after inflation. 
Alternatively, instead of being released from topological defects, X particles may have 
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been produced directly in the early Universe and, due to some unknown symmetries, have a 
very long lifetime comparable to the age of the Universe. In contrast to Weakly-Interacting 
Massive Particles (WIMPS) below a few hundred TeV which are the usual dark matter can- 
didates motivated by, for example, supersymmetry and can be produced by thermal freeze 
out, such superheavy X particles have to be produced non-thermally. Several such mech- 
anisms operating in the post- inflationary epoch in the early Universe have been studied. 
They include gravitational production through the effect of the expansion of the background 
metric on the vacuum quantum fluctuations of the X particle field, or creation during re- 
heating at the end of inflation if the X particle field couples to the inflaton field. The latter 
case can be divided into three subcases, namely "incoherent" production with an abundance 
proportional to the X particle annihilation cross section, non-adiabatic production in broad 
parametric resonances with the oscillating inflaton field during preheating (analogous to 
energy transfer in a system of coupled pendula), and creation in bubble wall collisions if 
inflation is completed by a first order phase transition. In all these cases, such particles, 
also called "WIMPZILLAs", would contribute to the dark matter and their decays could 
still contribute to UHE CR fluxes today, with an anisotropy pattern that reflects the dark 
matter distribution in the halo of our Galaxy. 

It is interesting to note that one of the prime motivations of the inflationary paradigm 
was to dilute excessive production of "dangerous relics" such as topological defects and su- 
perheavy stable particles. However, such objects can be produced right after inflation during 
reheating in cosmologically interesting abundances, and with a mass scale roughly given by 
the inflationary scale which in turn is fixed by the CMB anisotropies to ~ 10^'^ GeV [p8| . 
The reader will realize that this mass scale is somewhat above the highest energies observed 
in CRs, which implies that the decay products of these primordial relics could well have 
something to do with EHECRs which in turn can probe such scenarios! 

For dimensional reasons the spatially averaged X particle injection rate can only depend 
on the mass scale nix and on cosmic time t in the combination 



hx(t) = KUlxt' 



-4+p 



where k and p are dimensionless constants whose value depend on the specific top-down 
scenario |9^, For example, the case p = 1 is representative of scenarios involving release of 



X particles from topological defects, such as ordinary cosmic strings ||10CI|| , necklaces ||101| 



and magnetic monopoles ||102|| . This can be easily seen as follows: The energy density ps 
in a network of defects has to scale roughly as the critical density, ps oc pcrit oc where 
t is cosmic time, otherwise the defects would either start to overdose the Universe, or end 
up having a negligible contribution to the total energy density. In order to maintain this 
scaling, the defect network has to release energy with a rate given by ps = —aps/t oc t~^, 
where a = 1 in the radiation dominated aera, and a = 2/3 during matter domination. If most 
of this energy goes into emission of X particles, then typically k ~ (9(1). In the numerical 
simulations presented below, it was assumed that the X particles are nonrelativistic at decay. 

The X particles could be gauge bosons, Higgs bosons, superheavy fermions, etc. depending 
on the specific GUT. They would have a mass mx comparable to the symmetry breaking 
scale and would decay into leptons and/or quarks of roughly comparable energy. The quarks 
interact strongly and hadronize into nucleons (A^s) and pions, the latter decaying in turn 
into 7-rays, electrons, and neutrinos. Given the X particle production rate, dnx/dt, the 
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effective injection spectrum of particle species a {a = 'y,N,e^,i>) via the hadronic channel 
can be written as {dnx/dt){2/mx){dNa/dx), where x = 2E/mx-, and dNa/dx is the relevant 
fragmentation function (FF). 

We adopt the Local Parton Hadron Duality (LPHD) approximation ||103|] according to 
which the total hadronic FF, dNh/dx, is taken to be proportional to the spectrum of the 
partons (quarks/gluons) in the parton cascade (which is initiated by the quark through 
perturbative QCD processes) after evolving the parton cascade to a stage where the typical 
transverse momentum transfer in the QCD cascading processes has come down to ~ ~ 
few hundred MeV, where R is a. typical hadron size. The parton spectrum is obtained 
from solutions of the standard QCD evolution equations in modified leading logarithmic 



approximation (MLLA) which provides good fits to accelerator data at LEP energies ||103 



We will specifically use a recently suggested generalization of the MLLA spectrum that 
includes the effects of supersymmetry [ p.04|| . Within the LPHD hypothesis, the pions and 
nucleons after hadronization have essentially the same spectrum. The LPHD does not, 
however, fix the relative abundance of pions and nucleons after hadronization. Motivated 
by accelerator data, we assume the nucleon content of the hadrons to be in the range 3 
to 10%, and the rest pions distributed equally among the three charge states. According to 
recent Monte Carlo simulations ||1(J5||, the nucleon-to-pion ratio may be significantly higher 



in certain ranges of x values at the extremely high energies of interest here. Unfortunately, 
however, due to the very nature of these Monte Carlo calculations, it is difficult to understand 
the precise physical reason for the unexpectedly high baryon yield relative to mesons. While 
more of these Monte Carlo calculations of the relevant FFs in the future will hopefully clarify 
the situation, we will use here the range of /at ~ 3 to 10% mentioned above, which also seems 



to be supported by other recent Monte Carlo simulations ||106|| . The standard pion decay 



spectra then give the injection spectra of 7-rays, electrons, and neutrinos. For more details 
concerning uncertainties in the X particle decay spectra see Ref. [[107||. 



4.2 Numerical Simulations 

The 7-rays and electrons produced by X particle decay initiate electromagnetic (EM) cas- 
cades on low energy radiation fields such as the CMB. The high energy photons undergo 
electron-positron pair production (PP; 77b —>■ e~e~^), and at energies below ~ 10^"^ eV they 
interact mainly with the universal infrared and optical (IR/0) backgrounds, while above 
~ 100 EeV they interact mainly with the universal radio background (URB). In the Klein- 
Nishina regime, where the CM energy is large compared to the electron mass, one of the out- 
going particles usually carries most of the initial energy. This "leading" electron (positron) 
in turn can transfer almost all of its energy to a background photon via inverse Compton 
scattering (ICS; e% e'7). EM cascades are driven by this cycle of PP and ICS. The 
energy degradation of the "leading" particle in this cycle is slow, whereas the total number 
of particles grows exponentially with time. This makes a standard Monte Carlo treatment 
difficult. Implicit numerical schemes have therefore been used to solve the relevant kinetic 
equations. A detailed account of the transport equation approach used in the calculations 



whose results are presented in this contribution can be found in Ref. ||108|| . All EM interac- 
tions that influence the 7-ray spectrum in the energy range 10^ eV < E < 10^^ eV, namely 
PP, ICS, triplet pair production (TPP; 6% ee~e~^), and double pair production (DPP, 
776 e~e^e~e^), as well as synchrotron losses of electrons in the large scale extragalactic 
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magnetic field (EGMF), are included. 

Similarly to photons, UHE neutrinos give rise to neutrino cascades in the primordial 
neutrino background via exchange of W and Z bosons [|109| , |11CI|] . Besides the secondary 
neutrinos which drive the neutrino cascade, the W and Z decay products include charged 
leptons and quarks which in turn feed into the EM and hadronic channels. Neutrino inter- 
actions become especially significant if the relic neutrinos have masses mj, in the eV range 
and thus constitute hot dark matter, because the Z boson resonance then occurs at an UHE 
neutrino energy Ej-cs = 4 x 10'^^{eV /m,^) eV. In fact, this has been proposed as a significant 



source of EHECRs ||111| , |112|| . Motivated by recent experimental evidence for neutrino mass 
we assumed a mass of 1 eV for all three neutrino flavors (for simplicity) and implemented 
the relevant W boson interactions in the t-channel and the Z boson exchange via t- and 
s-channel. Hot dark matter is also expected to cluster, potentially increasing secondary 
7-ray and nucleon production ||111| , |112|| . This influences mostly scenarios where X decays 
into neutrinos only. We parametrize massive neutrino clustering by a length scale li, and an 
overdensity over the average density ri^. The Fermi distribution with a velocity dispersion 
V yields f,y < v^ml /(27r)3/2/n^ ~ 330 (i;/500km sec {rriiy / eV)^ ||114|| . Therefore, values of 



~ few Mpc and fi, ~ 20 are conceivable on the local Supercluster scale ||112 |. 

The relevant nucleon interactions implemented are pair production by protons {pyi, 
pe~e~^), photoproduction of single or multiple pions {N% N rni, n > 1), and neutron 
decay. In TD scenarios, the particle injection spectrum is generally dominated by the "pri- 
mary" 7-rays and neutrinos over nucleons. These primary 7-rays and neutrinos are produced 
by the decay of the primary pions resulting from the hadronization of quarks that come from 
the decay of the X particles. The contribution of secondary 7-rays, electrons, and neutrinos 
from decaying pions that are subsequently produced by the interactions of nucleons with 
the CMB, is in general negligible compared to that of the primary particles; we nevertheless 
include the contribution of the secondary particles in our code. 

In principle, new interactions such as the ones involving a TeV quantum gravity scale 
can not only modify the interactions of primary particles in the detector, as discussed in 
Sect. 3.1, but also their propagation. However, 7— rays and nuclei interact mostly with 
the CMB and IR for which the CM energy is at most ~ 30(^/10^^ GeV)^/^ Q^y ^t such 
energies the new interactions are much weaker than the dominating electromagnetic and 
strong interactions. It has been suggested recently [|113|| that new interactions may notably 
influence UHE neutrino propagation. However, for neutrinos of mass m^, the CM energy is 
~ 100(i?/10^^ GeV)^/2(m^/0.1 eV)^/2 GeV, and therefore, UHE neutrino propagation would 
only be significantly modified for neutrino masses significantly larger than 0.1 eV. We will 
therefore ignore this possibility here. 

We assume a flat Universe with no cosmological constant, and a Hubble constant of 
h = 0.65 in units of 100 km sec~^Mpc~^ throughout. The numerical calculations follow 
all produced particles in the EM, hadronic, and neutrino channel, whereas the often-used 
continuous energy loss (CEL) approximation (e.g., ||115|| ) follows only the leading cascade 
particles. The CEL approximation can significantly underestimate the cascade flux at lower 
energies. 

The two major uncertainties in the particle transport are the intensity and spectrum of 
the URB for which there exists only an estimate above a few MHz frequency ||116|| , and the 
average value of the EGMF. To bracket these uncertainties, simulations have been performed 
for the observational URB estimate from Ref. |p.l6|| that has a low-frequency cutoff at 2 MHz 
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("minimal"), and the medium and maximal theoretical estimates from Ref. ||117]| , as well as 
for EGMFs between zero and 10~^ G, the latter motivated by limits from Faraday rotation 
measurements, see Sect. 5.2 below. A strong URB tends to suppress the UHE 7-ray flux by 
direct absorption whereas a strong EGMF blocks EM cascading (which otherwise develops 
efficiently especially in a low URB) by synchrotron cooling of the electrons. For the IR/0 
background we used the most recent data [|118||. 



4.3 Results: 7— ray and Nucleon Fluxes 
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Figure 4: Predictions for the differential fluxes of 7— rays (solid line) and protons and neutrons (dotted 
line) in a TD model characterized by p = 1, rax = 10^^ GeV, and the decay mode X ^ q + q, assuming 



the supersymmetric modification of the fragmentation function |104|, with a fraction of about 10% nucleons. 
The defects have been assumed to be homogeneously distributed. The calculation used the code described in 
Ref. [ [107| and assumed the strongest URB version from Ref. [117| and an EGMF <C 10~^^ G. 1 sigma error 
bars are the combined data from the Haverah Park Q, the Fly's Eye and the AGASA |^ experiments 
above 10^^ eV. Also shown are piecewise power law fits to the observed charged GR flux (thick solid line) and 
the EGRET measurement of the diffuse 7— ray flux between 30 MeV and 100 GeV [ |119| (solid line on left 
margin). Points with arrows represent upper limits on the 7— ray flux from the HEGRA, the Utah-Michigan, 
the EAS-TOP, and the CASA-MIA experiments, as indicated. 



Fig. ^ shows results from Ref. | ]107| ] for the time averaged 7— ray and nucleon fluxes in 
a typical TD scenario, assuming no EGMF, along with current observational constraints on 
the 7— ray flux. The spectrum was optimally normalized to allow for an explanation of the 
observed EHECR events, assuming their consistency with a nucleon or 7— ray primary. The 
flux below < 2 X 10^^ eV is presumably due to conventional acceleration in astrophysical 
sources and was not flt. Similar spectral shapes have been obtained in Ref. [|120|] , where 
the normalization was chosen to match the observed differential flux at 3 x 10^*^ eV. This 
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Figure 5: Same as Fig. |, but for an EGMF of lO'^ G. 



normalization, however, leads to an overproduction of the integral flux at higher energies, 
whereas above 10^° eV, the fits shown in Figs. ^ and ^ have likelihood significances above 
50% (see Ref. ||121| | for details) and are consistent with the integral flux above 3 x 10^° eV 
estimated in Refs. [|, ^. The PP process on the CMB depletes the photon flux above 100 
TeV, and the same process on the IR/0 background causes depletion of the photon flux 
in the range 100 GeV-100 TeV, recycling the absorbed energies to energies below 100 GeV 
through EM cascading (see Fig. 1). The predicted background is not very sensitive to the 
specific IR/0 background model, however ||122|| . The scenario in Fig. ^obviously obeys all 
current constraints within the normalization ambiguities and is therefore quite viable. Note 
that the diffuse 7— ray background measured by EGRET ||119|| up to 10 GeV puts a strong 
constraint on these scenarios, especially if there is already a significant contribution to this 
background from conventional sources such as unresolved 7— ray blazars [123]. However, the 
7— ray background constraint can be circumvented by assuming that TDs or the decaying 
long lived X particles do not have a uniform density throughout the Universe but cluster 
within galaxies ||124||. As can also be seen, at energies above 100 GeV, TD models are not 



significantly constrained by observed 7— ray fluxes yet (see Ref. |]T5| for more details on these 
measurements) . 

Fig. 1^ shows results for the same TD scenario as in Fig. ^, but for a high EGMF ~ 10~^ G, 
somewhat below the current upper limit, see Eq. (|l^) below. In this case, rapid synchrotron 
cooling of the initial cascade pairs quickly transfers energy out of the UHE range. The UHE 
7— ray flux then depends mainly on the absorption length due to pair production and is 
typically much lower ||115| , p.25|| . (Note, though, that for mx ^ 10^^ eV, the synchrotron 
radiation from these pairs can be above lO^'' eV, and the UHE flux is then not as low as one 
might expect.) We note, however, that the constraints from the EGRET measurements do 
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not change significantly with the EGMF strength as long as the nucleon flux is comparable 
to the 7— ray flux at the highest energies, as is the case in Figs. | and ^. The results of 



Ref. 1 107 1 differ from those of Ref. [120| which obtained more stringent constraints on TD 



models because of the use of an older fragmentation function from Ref. ||126|| , and a stronger 
dependence on the EGMF because of the use of a weaker EGMF which lead to a dominance 
of 7— rays above ~ lO^'' eV. 

The energy loss and absorption lengths for UHE nucleons and photons are short (< 100 
Mpc). Thus, their predicted UHE fluxes are independent of cosmological evolution. The 
7— ray flux below ~ 10^^ eV, however, scales as the total X particle energy release integrated 
over all redshifts and increases with decreasing p ||127|| . For mx = 2 x 10^^ GeV, scenarios 
with p < 1 are therefore ruled out (as can be inferred from Figs. ^ and P), whereas constant 
comoving injection models {p = 2) are well within the limits. 

We now turn to signatures of TD models at UHE. The full cascade calculations predict 
7— ray fluxes below 100 EeV that are a factor ~ 3 and ~ 10 higher than those obtained 
using the GEL or absorption approximation often used in the literature, in the case of strong 
and weak URB, respectively. Again, this shows the importance of non- leading particles in 
the development of unsaturated EM cascades at energies below ~ 10^^ eV. Our numerical 
simulations give a 7/GR flux ratio at 10^^ eV of ~ 0.1. The experimental exposure required 
to detect a 7— ray flux at that level is ~ 4 x 10^^ cm^ sec sr, about a factor 10 smaller than 
the current total experimental exposure. These exposures are well within reach of the Pierre 
Auger Gosmic Ray Observatories , which may be able to detect a neutral GR component 
down to a level of 1% of the total flux. In contrast, if the EGMF exceeds ~ 10^^^ G, then 
UHE cascading is inhibited, resulting in a lower UHE 7— ray spectrum. In the 10^^ G 
scenario of Fig. ^, the 7/GR flux ratio at 10^^ eV is 0.02, signiflcantly lower than for no 
EGMF. 

It is clear from the above discussions that the predicted particle fluxes in the TD scenario 
are currently uncertain to a large extent due to particle physics uncertainties (e.g., mass and 
decay modes of the X particles, the quark fragmentation function, the nucleon fraction /^r, 
and so on) as well as astrophysical uncertainties (e.g., strengths of the radio and infrared 
backgrounds, extragalactic magnetic flelds, etc.). More details on the dependence of the 
predicted UHE particle spectra and composition on these particle physics and astrophysical 



uncertainties are contained in Ref. [|107|| . A detailed study of the uncertainties involved in 
the propagation of UHE nucleons, 7— rays, and neutrinos is currently underway [|128|] . 

We stress here that there are viable TD scenarios which predict nucleon fluxes that 
are comparable to or even higher than the 7— ray flux at all energies, even though 7— rays 
dominate at production. This occurs, e.g., in the case of high URB and/or for a strong 
EGMF, and a nucleon fragmentation fraction of ~ 10%; see, for example. Fig. ^ Some 
of these TD scenarios would therefore remain viable even if EHEGR induced EAS should 
be proven inconsistent with photon primaries (see, e.g., Ref. ||129|| ). This is in contrast 
to scenarios with decaying massive dark matter in the Galactic halo which, due to the 
lack of absorption, predict compositions directly given by the fragmentation function, i.e. 
domination by 7— rays. 

The normalization procedure to the EHEGR flux described above imposes the constraint 
Qehecr ~ 10~^^ eV cm~'^ sec"^ within a factor of a few [|120| , [lOT] , |130|| for the total energy 



release rate Qq from TDs at the current epoch. In most TD models, because of the unknown 
values of the parameters involved, it is currently not possible to calculate the exact value 
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of Qo from first principles, although it has been shown that the required values of Qq (in 
order to explain the EHECR flux) mentioned above are quite possible for certain kinds of 
TDs. Some cosmic string simulations and the necklace scenario suggest that defects may 
lose most of their energy in the form of X particles and estimates of this rate have been 
given [|131| , |101|| . If that is the case, the constraint on Qehecr translates via Eq. (H) into 
a limit on the symmetry breaking scale r] and hence on the mass mx of the X particle: 
1] ~ mx ^ 10^^ GeV ||132|| . Independently of whether or not this scenario explains EHECR, 



the EGRET measurement of the diffuse GeV 7— ray background leads to a similar bound, 
Qem ~ 2.2 X 10^^'^ /;,(3p — 1) eVcm"'^ sec~^, which leaves the bound on rj and mx practically 
unchanged. Furthermore, constraints from limits on CMB distortions and light element 
abundances from ^He-photodisintegration are comparable to the bound from the directly 
observed diffuse GeV 7-rays | 127| |. That these crude normalizations lead to values of r/ in the 



right range suggests that defect models require less fine tuning than decay rates in scenarios 
of metastable massive dark matter. 



4.4 Results: Neutrino Fluxes 

As discussed in Sect. 4.1, in TD scenarios most of the energy is released in the form of 
EM particles and neutrinos. If the X particles decay into a quark and a lepton, the quark 
hadronizes mostly into pions and the ratio of energy release into the neutrino versus EM 
channel is r ~ 0.3. 

Fig. ^ shows predictions of the total neutrino flux for the same TD model on which Fig. ^ 
is based. In the absence of neutrino oscillations the electron neutrino and anti-neutrino 
fluxes are about a factor of 2 smaller than the muon neutrino and anti-neutrino fluxes, 
whereas the r— neutrino flux is in general negligible. In contrast, if the interpretation of the 
atmospheric neutrino deficit in terms of nearly maximal mixing of muon and r— neutrinos 
proves correct, the muon neutrino fluxes shown in Fig. ^ would be maximally mixed with 
the r— neutrino fluxes. To put the TD component of the neutrino flux in perspective with 



contributions from other sources. Fig. M also shows the atmospheric neutrino flux |135 



typical prediction for the diffuse flux from photon optically thick proton blazars |136| that 
are not subject to the Waxman Bahcall bound and were normalized to recent estimates of 
the blazar contribution to the diffuse 7— ray background [p.23|| , and the flux range expected 
for "cosmogenic" neutrinos created as secondaries from the decay of charged pions produced 
by UHE nucleons ||137||. The TD flux component clearly dominates above ~ 10^^ eV. 



In order to translate neutrino fluxes into event rates, one has to fold in the interaction 
cross sections with matter. At UHEs these cross sections are not directly accessible to 
laboratory measurements. Resulting uncertainties therefore translate directly to bounds on 
neutrino fluxes derived from, for example, the non-detection of UHE muons produced in 
charged- current interactions. In the following, we will assume the estimate Eq. (|^) based on 
the Standard Model for the charged-current muon-neutrino-nucleon cross section ai,x if not 
indicated otherwise. 

For an (energy dependent) ice or water equivalent acceptance A{E) (in units of volume 
times solid angle), one can obtain an approximate expected rate of UHE muons produced 
by neutrinos with energy > E, R{E), by multiplying A{E)ai,j\f{E)nu20 (where nHjO is the 
nucleon density in water) with the integral muon neutrino flux ~ Ej^,^. This can be used to 
derive upper limits on diffuse neutrino fluxes from a non-detection of muon induced events. 
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Figure 6: Neutrino fluxes predicted by the TD model from fig. ^, marked by "SLB Y98" [ jl07[ . Furthermore 
shown are experimental neutrino flux limits from the Frejus underground detector [133|, the Fly's Eye JtO] , 
the Goldstone radio telescope [134|, and the Antarctic Muon and Neutrino Detector Array (AMANDA) 
neutrino telescope as well as projected neutrino flux sensitivities of ICECUBE, the planned kilometer 
scale extension of AMANDA the Pierre Auger Project |Q (for electron and tau neutrinos separately) 
and the proposed space based OWL [^2| concept. For comparison also shown are the atmospheric neutrino 
background |l35| (hatched region marked "atmospheric"), and neutrino flux predictions for a model of AGN 
optically thick to nucleons ("AGN"), and for UHECR interactions with the CMB (^S^l ("iV7", dashed 
range indicating typical uncertainties for moderate source evolution). The top-down fluxes are shown for 
electron-, muon, and tau- neutrinos separately, assuming no (lower zy,--curve) and maximal — Vr mixing 
(upper i/T--curve, which would then equal the ^-^-flux), respectively. 



Fig. ^ shows bounds obtained from several experiments: The Frejus experiment derived upper 
bounds for E > 10^^ eV from their non-detection of almost horizontal muons with an energy 
loss inside the detector of more than 140 MeV per radiation length ||133|] . The EAS-TOP 
collaboration published two limits from horizontal showers, one in the regime 10^^ — 10^^ eV, 
where non-resonant neutrino-nucleon processes dominate, and one at the Glashow resonance 
which actually only applies to Ue ||138|| . The Fly's Eye experiment derived upper bounds for 

1020 eV 



70 1 from the non-observation of deeply 



the energy range between ~ 10 eV and 
penetrating particles. The AKENO group has published an upper bound on the rate of near- 
horizontal, muon-poor air showers |jl39|| . Horizontal air showers created by electrons, muons 
or tau leptons that are in turn produced by charged-current reactions of electron, muon or 
tau neutrinos within the atmosphere have recently also been pointed out as an important 



method to constrain or measure UHE neutrino fluxes [53| with next generation detectors. 



The p = TD model BHSO from the early work of Ref. ||99| is not only ruled out by the 
constraints from Sect. 4.3, but also by some of the experimental limits on the UHE neutrino 
flux, as can be seen in Fig. ^. Further, although both the BHSl and the SLBY98 models 
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correspond to p = 1, the UHE neutrino flux above ~ 10^° eV in the latter is almost two 
orders of magnitude smaller than in the former. The main reason for this is the different flux 
normalization adopted in the two papers: First, the BHSl model was obtained by normalizing 
the predicted proton flux to the observed UHECR flux at ~ 4 x 10^^ eV, whereas in the 
SLBY98 model the actually "visible" sum of the nucleon and 7— ray fluxes was normalized in 



an optimal way. Second, the BHSl assumed a nucleon fraction about a factor 3 smaller [99 



Third, the BHSl scenario used an older fragmentation function from Ref. ||126|| which has 
more power at larger energies. Clearly, the SLBY98 model is not only consistent with the 
constraints discussed in Sect. 4.3, but also with all existing neutrino flux limits within 2-3 
orders of magnitude. 

What, then, are the prospects of detecting UHE neutrino fluxes predicted by TD models? 
In a 1 km'^ 27r sr size detector, the SLBY98 scenario from Fig. ^, for example, predicts a muon- 
neutrino event rate of ~ 0.15 yr^^, and an electron neutrino event rate of ~ 0.089 yr^^ above 
10^^ eV, where "backgrounds" from conventional sources should be negligible. Further, the 
muon-neutrino event rate above 1 PeV should be ~ 1.2 yr~^, which could be interesting 
if conventional sources produce neutrinos at a much smaller flux level. Of course, above ~ 
100 TeV, instruments using ice or water as detector medium, have to look at downward going 
muon and electron events due to neutrino absorption in the Earth. However, r— neutrinos 



obliterate this Earth shadowing effect due to their regeneration from r decays ||140|| . The 



presence of r— neutrinos, for example, due to mixing with muon neutrinos, as suggested 
by recent experimental results from Super-Kamiokande, can therefore lead to an increased 
upward going event rate ||141||. For recent compilations of UHE neutrino flux predictions 



from astrophysical and TD sources see Refs. |142|| and references therein. 

For detectors based on the fluorescence technique such as the HiRes and the Telescope 
Array |]39| (see Sect. 2), the sensitivity to UHE neutrinos is often expressed in terms of 
an effective aperture a{E) which is related to A{E) by a{E) = A{E)a^N{E)nii^o. For 
the cross section of Eq. (^, the apertures given in Ref. for the HiRes correspond to 
A(E) ~ 3km^ X 27rsr for E > 10^^ eV for muon neutrinos. The expected acceptance of the 
ground array component of the Pierre Auger project for horizontal UHE neutrino induced 
events is y4(10^^eV) ~ 20km^sr and y4(10^^eV) ~ 200km^sr with a duty cycle close 
to 100%. We conclude that detection of neutrino fluxes predicted by scenarios such as the 
SLBY98 scenario shown in Fig. ^ requires running a detector of acceptance > 10 km^ x 2n sr 
over a period of a few years. Apart from optical detection in air, water, or ice, other methods 
such as acoustical and radio detection pSf (see, e.g., the RICE project |^ for the latter) or 



even detection from space |4|, ^ appear to be interesting possibilities for detection 

concepts operating at such scales (see Sect. 2). For example, the space based OWL/AirWatch 
satellite concept would have an aperture of ~ 3 x 10^ km^ sr in the atmosphere, corresponding 
to A{E) ~ 6 X lO^km^sr for E > lO^^eV, with a duty cycle of ~ 0.08 H 11. The 
backgrounds seem to be in general negligible ||110| , |143|| . As indicated by the numbers above 
and by the projected sensitivities shown in Fig. ^, the Pierre Auger Project and especially 
the space based AirWatch type projects should be capable of detecting typical TD neutrino 
fluxes. This applies to any detector of acceptance > lOOkm^sr. Furthermore, a 100 day 
search with a radio telescope of the NASA Goldstone type for pulsed radio emission from 
cascades induced by neutrinos or cosmic rays in the lunar regolith could reach a sensitivity 
comparable or better to the Pierre Auger sensitivity above ~ 10^^ eV 



A more model independent estimate ||130|| for the average event rate R{E) can be made 
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if the underlying scenario is consistent with observational nucleon and 7— ray fluxes and the 
bulk of the energy is released above the PP threshold on the CMB. Let us assume that the 
ratio of energy injected into the neutrino versus EM channel is a constant r. As discussed in 
Sect. 4.3, cascading effectively reprocesses most of the injected EM energy into low energy 
photons whose spectrum peaks at ~ 10 GeV ||122|| . Since the ratio r remains roughly un- 
changed during propagation, the height of the corresponding peak in the neutrino spectrum 
should roughly be r times the height of the low-energy 7— ray peak, i.e., we have the condition 
max^; E^j^^{E) ~ r max^ [E'^j^{E)] . Imposing the observational upper limit on the diffuse 



7— ray flux around 10 GeV shown in Fig. ^, max^; E'^j^, {E) 



< 2 X lO^reVcm ^sec ^sr" 



then bounds the average diffuse neutrino rate above PP threshold on the CMB, giving 



R{E) < 0.34 r 



A{E) 



1 km X 2tt ST 



E 



1019 eV 



-0.6 



yr-i {E>10^^eV) 



(9) 



For r < 20{E/10^^ eV)^-^ this bound is consistent with the flux bounds shown in Fig. ^that 
are dominated by the Fly's Eye constraint at UHE. We stress again that TD models are not 
subject to the Waxman Bahcall bound because the nucleons produced are considerably less 
abundant than and are not the primaries of produced 7— rays and neutrinos. 




E (eV) 

Figure 7: Flux predictions for a TD model characterized by p = 1, mx = 10^'* GeV, with X particles 
exclusively decaying into neutrino-antineutrino pairs of all flavors (with equal branching ratio), assuming a 
neutrino mass mi, = 1 eV. For neutrino clustering, an overdensity of ~ 3 ove r a scale of ~ 5 Mpc was 
assumed. The calculation assumed the strongest URB version from Ref. |ll7| and an EGMF ^ lO^^^G. 
The line key is as in Figs. ^ and ^. 

In typical TD models such as the one discussed above where primary neutrinos are pro- 
duced by pion decay, r ^ 0.3. However, in TD scenarios with r ^ 1 neutrino fluxes are only 
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limited by the condition that the secondary 7— ray flux produced by neutrino interactions 
with the relic neutrino background be below the experimental limits. An example for such a 
scenario is given by X particles exclusively decaying into neutrinos (although this is not very 
likely in most particle physics models, but see Ref. ||107|| and Fig. |^ for a scenario involving 
topological defects and Ref. ||145|| for a scenario involving decaying superheavy relic particles, 
both of which explain the observed EHECR events as secondaries of neutrinos interacting 
with the primordial neutrino background). Such scenarios predict appreciable event rates 
above ~ 10^^ eV in a km^ scale detector, but require unrealistically strong clustering of relic 
neutrinos (a homogeneous relic neutrino overdensity would make the EGRET constraint 
only more severe because neutrino interactions beyond ~ 50 Mpc contribute to the GeV 
7— ray background but not to the UHECR flux). A detection would thus open the exciting 
possibility to establish an experimental lower limit on r. Being based solely on energy con- 
servation, Eq. (|^) holds regardless of whether or not the underlying TD mechanism explains 
the observed EHECR events. 

The transient neutrino event rate could be much higher than Eq. (|^) in the direction to 
discrete sources which emit particles in bursts. Corresponding pulses in the EHE nucleon 
and 7— ray fluxes would only occur for sources nearer than ~ 100 Mpc and, in case of 
protons, would be delayed and dispersed by deflection in Galactic and extragalactic magnetic 
fields ||146| , |147|| . The recent observation of a possible clustering of CRs above ~ 4 x 10^^ eV by 
the AGASA experiment [|148|| might suggest sources which burst on a time scale tb 1 yr. A 
burst fluence of ~ r A{E) /Ikm!^ x 27rsr eV)~^'^ neutrino induced events within a 

time tf, could then be expected. Associated pulses could also be observable in the GeV — TeV 
7— ray flux if the EGMF is smaller than ~ 10"^^ G in a signiflcant fraction of extragalactic 
space 
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In contrast to roughly homogeneous sources and/or mechanisms with branching ratios 
r 3> 1, in scenarios involving clustered sources such as metastable superheavy relic particles 
decaying with r ~ 1, the neutrino flux is comparable to (not signiflcantly larger than) 
the UHE photon plus nucleon fluxes. This can be understood because the neutrino flux 
is dominated by the extragalactic contribution which scales with the extragalactic nucleon 
and 7— ray contribution in exactly the same way as in the unclustered case, whereas the 
extragalactic contribution to the "visible" flux to be normalized to the UHECR data is 
much smaller in the clustered case. The resulting neutrino fluxes would be hardly detectable 
even with next generation experiments. 



5 UHE Cosmic Rays and Cosmological Large Scale 
Magnetic Fields 

5.1 Deflection and Delay of Charged Hadrons 

Whereas for UHE electrons the dominant influence of large scale magnetic flelds is syn- 
chrotron loss rather than deflection, for charged hadrons the opposite is the case. A rela- 
tivistic particle of charge qe and energy E has a gyroradius Vg ~ E/{qeB±) where 5^ is 
the fleld component perpendicular to the particle momentum. If this fleld is constant over a 
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distance d, this leads to a deflection angle 

i^.0,5r,(^)-'(^)(^). ,10) 

Magnetic fields beyond the Galactic disk are poorly known and include a possible ex- 
tended field in the halo of our Galaxy and a large scale EGMF. In both cases, the magnetic 
field is often characterized by an r.m.s. strength B and a correlation length Ic, i.e. it is 
assumed that its power spectrum has a cut-off in wavenumber space at = 27r/Zc and in real 
space it is smooth on scales below Ic- If we neglect energy loss processes for the moment, 
then the r.m.s. deflection angle over a distance d in such a field is 9{E,d) ~ {2dlc/9Y^'^ /rg, 
or 

OiE, ,1)^0.8' J ( ^] ( ^) ( J^) . (11) 



1020 eV/ VlOMpc/ VlMpc/ V10-9G 

for d > Ic, where the numerical prefactors were calculated from the analytical treatment in 
Ref. ||146|| . There it was also pointed out that there are two different limits to distinguish: 
For d9{E, d) ^ Ic, particles of all energies "see" the same magnetic field realization during 
their propagation from a discrete source to the observer. In this case, Eq. (|1T]) gives the 
typical coherent deflection from the line-of-sight source direction, and the spread in arrival 
directions of particles of different energies is much smaller. In contrast, for d6{E, d) ^ Ic, 
the image of the source is washed out over a typical angular extent again given by Eq. ([TT|), 
but in this case it is centered on the true source direction. If d9{E, d) ~ Ic, the source may 
even have several images, similar to the case of gravitational lensing. Therefore, observing 
images of UHECR sources and identifying counterparts in other wavelengths would allow 
one to distinguish these limits and thus obtain information on cosmic magnetic flelds. If d is 
comparable to or larger than the interaction length for stochastic energy loss due to photo- 
pion production or photodisintegration, the spread in deflection angles is always comparable 
to the average deflection angle. 

Deflection also implies an average time delay of t{E, d) ^ d6{E, d^/A, or 

3 2 f \ ^ f ^ \ f \ f B 



r(B.d).1.5xlO^,^^_j _ U_) (12) 



relative to rectilinear propagation with the speed of light. It was pointed out in Ref. ||150 
that, as a consequence, the observed UHECR spectrum of a bursting source at a given 
time can be different from its long-time average and would typically peak around an energy 
Eq, given by equating t{E, d) with the time of observation relative to the time of arrival 
for vanishing time delay. Higher energy particles would have passed the observer already, 
whereas lower energy particles would not have arrived yet. Similarly to the behavior of 
deflection angles, the width of the spectrum around Eq would be much smaller than £"0 if 
both d is smaller than the interaction length for stochastic energy loss and d6{E, d) -C Ic- In 
all other cases the width would be comparable to Eq- 

Constraints on magnetic fields from defiection and time delay cannot be studied sepa- 
rately from the characteristics of the "probes" , namely the UHECR sources, at least as long 
as their nature is unknown. An approach to the general case is discussed in Sect. 5.3. 
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5.2 Constraints on EHECR Source Locations 

As pointed out in Sect. 1, nucleons, nuclei, and 7— rays above a few 10^^ eV cannot have 
originated much further away than ^ 50 Mpc. Together with Eq. (0) this imphes that 
above a few 10^^ eV the arrival direction of such particles should in general point back to their 
source within a few degrees |T7]. This argument is often made in the literature and follows 



from the Faraday rotation bound on the EGMF and a possible extended field in the halo of 
our Galaxy, which in its historical form reads Bl^"^ < 10~^GMpc^''^ ||151|| , as well as from 



the known strength and scale height of the field in the disk of our Galaxy, i?^ ~ 3 x 10~^ G, 
Ig < 1 kpc. Furthermore, the deflection in the disk of our Galaxy can be corrected for 
in order to reconstruct the extragalactic arrival direction: Maps of such corrections as a 
function of arrival direction have been calculated in Refs. ||152|| for plausible models of the 
Galactic magnetic field. The deflection of UHECR trajectories in the Galactic magnetic field 
may, however, also give rise to several other important effects ||153[ such as (de)magnification 



of the UHECR fluxes due to the magnetic lensing effect mentioned in the previous section 
(which can modify the UHECR spectrum from individual sources), formation of multiple 
images of a source, and apparent "blindness" of the Earth towards certain regions of the sky 
with regard to UHECRs. These effects may in turn have important implications for UHECR 



source locations. In fact, it was recently claimed ||154|] that, assuming a certain model of 
the magnetic fields in the galactic winds, the highest energy cosmic ray events could all 
have originated in the Virgo cluster or specifically in the radio galaxy M87. However, as 



was subsequently pointed out in Ref. ||155|| , this galactic wind model leads to focusing of all 



positively charged highest energy particles to the North galactic pole and, consequently, this 
can not be interpreted as evidence for a point source situated close to the North galactic 
pole. 

However, important modifications of the Faraday rotation bound on the EGMF have 
recently been discussed in the literature: The average electron density which enters estimates 
of the EGMF from rotation measures, can now be more reliably estimated from the baryon 
density flbh^ ~ 0.02, whereas in the original bound the closure density was used. Assuming 
an unstructured Universe and fig = 1 results in the much weaker bound ||156| 



5<3xlO-M-^ \\^—\ G, 13 

VO-02 / VO-65/ VMpc; ' ^ ' 

which suggests much stronger deflection. However, taking into account the large scale struc- 
ture of the Universe in the form of voids, sheets, filaments etc., and assuming flux freezing 
of the magnetic fields whose strength then approximately scales with the 2/3 power of the 
local density, leads to more stringent bounds: Using the Lyman a forest to model the density 
distribution yields ||156 

5<10"^-10"^G (14) 

for the large scale EGMF for coherence scales between the Hubble scale and 1 Mpc. This 
estimate is closer to the original Faraday rotation limit. However, in this scenario the 



maximal fields in the sheets and voids can be as high G I153, |T56| . 

Therefore, according to Eq. ([TI|) and (|14D, deflection of UHECR nucleons is still ex- 
pected to be on the degree scale if the local large scale structure around the Earth is not 
strongly magnetized. However, rather strong deflection can occur if the Supergalactic Plane 
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is strongly magnetized, for particles originating in nearby galaxy clusters where magnetic 
fields can be as high as 10~^ G |p.51|| (see Sect. 5.3 below) and/or for heavy nuclei such as 
In this case, magnetic lensing in the EGMF can also play an important role in 



iron 
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determining UHECR source locations [15S, 159 



5.3 Angle-Energy-Time Images of UHECR Sources 
Small Deflection 

For small deflection angles and if photo-pion production is important, one has to resort to 
numerical Monte Carlo simulations in 3 dimensions. Such simulations have been performed 



in Ref. ||160|| for the case d6{E, d) ^ Ic and in Refs. |p.47|, |161|, |162|| for the general case 



T 




T [yr-] 



Fi gure 8: Contour plot of the UHECR image of a bursting sour ce at d — 30Mpc, projected onto the 
time-energy plane, with B = 2 x 10^^" G, Ic — I Mpc, from Ref. |147|. The contours decrease in steps of 0.2 
in the logarithm to base 10. The dotted line indicates the energy-time delay correlation t{E, d) oc E~'^ as 
would be obtained in the absence of pion production losses. Clearly, d9{E, d) ^ Ic in this example, since for 
E < 4 X 10^^ eV, the width of the energy distribution at any given time is much smaller than the average 
(see Sect. 5.1). The dashed lines, which are not resolved here, indicate the location (arbitrarily chosen) of 
the observational window, of length Tohs = 5 yr. 



In Refs. ||147| , |161| , |162|| the Monte Carlo simulations were performed in the following 
way: The magnetic field was represented as a Gaussian random field with zero mean and a 
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power spectrum with {B'^{k)) oc k""" for k < and {B'^{k)) = otherwise, where k^ = 27r//c 
characterizes the numerical cut-off scale and the r.m.s. strength is = /q°° dkk^ {B'^{k)). 
The field is then calculated on a grid in real space via Fourier transformation. For a given 
magnetic field realization and source, nucleons with a uniform logarithmic distribution of 
injection energies are propagated between two given points (source and observer) on the grid. 
This is done by solving the equations of motion in the magnetic field interpolated between 
the grid points, and subjecting nucleons to stochastic production of pions and (in case of 
protons) continuous loss of energy due to PP. Upon arrival, injection and detection energy, 
and time and direction of arrival are recorded. From many (typically 40000) propagated 
particles, a histogram of average number of particles detected as a function of time and 
energy of arrival is constructed for any given injection spectrum by weighting the injection 
energies correspondingly. This histogram can be scaled to any desired total fiuence at the 
detector and, by convolution in time, can be constructed for arbitrary emission time scales 
of the source. An example for the distribution of arrival times and energies of UHECRs from 
a bursting source is given in Fig. ^ 

We adopt the following notation for the parameters: rioo denotes the time delay due to 
magnetic deflection ai E = 100 EeV and is given by Eq. (|T2|) in terms of the magnetic field 
parameters; T5 denotes the emission time scale of the source; Ts ^ lyr corresponds to a 
burst, and Tg ^ lyr (roughly speaking) to a continuous source; 7 is the differential index 
of the injection energy spectrum; Nq denotes the fiuence of the source with respect to the 
detector, i.e., the total number of particles that the detector would detect from the source 
on an infinite time scale; finally, C is the likelihood function of the above parameters. 

By putting windows of width equal to the time scale of observation over these histograms 
one obtains expected distributions of events in energy and time and direction of arrival for 
a given magnetic field realization, source distance and position, emission time scale, total 
fiuence, and injection spectrum. Examples of the resulting energy spectrum are shown in 
Fig. 1^. By dialing Poisson statistics on such distributions, one can simulate corresponding 
observable event clusters. 

Conversely, for any given real or simulated event cluster, one can construct a likelihood 
of the observation as a function of the time delay, the emission time scale, the differential 
injection spectrum index, the fiuence, and the distance. In order to do so, and to obtain the 
maximum of the likelihood, one constructs histograms for many different parameter com- 
binations as described above, randomly puts observing time windows over the histograms, 
calculates the likelihood function from the part of the histogram within the window and the 
cluster events, and averages over different window locations and magnetic field realizations. 

In Ref. ||161|| this approach has been applied to and discussed in detail for the three pairs 



observed by the AGASA experiment ||148|| , under the assumption that all events within a 



pair were produced by the same discrete source. Although the inferred angle between the 
momenta of the paired events acquired in the EGMF is several degrees |p.6^], this is not 



necessarily evidence against a common source, given the uncertainties in the Galactic field 
and the angular resolution of AGASA which is ~ 2.5°. As a result of the likelihood analysis, 
these pairs do not seem to follow a common characteristic; one of them seems to favor a 
burst, another one seems to be more consistent with a continuously emitting source. The 
current data, therefore, does not allow one to rule out any of the models of UHECR sources. 
Furthermore, two of the three pairs are insensitive to the time delay. However, the pair 
which contains the 200 EeV event seems to significantly favor a comparatively small average 
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Figure 9: Energy spectra for a continuous source (solid line), and for a burst (dashed line), from Ref. |147[ . 
Both spectra are normalized to a total of 50 particles detected. The parameters corresponding to the 
continuous source case are: Ts — 10^ yr, noo = 1-3 x 10'^ yr, and the time of observation is i = 9 x 10"^ yr, 
relative to rectilinear propagation with the speed of light. A low energy cutoff results at the energy Es = 
4 X 10^^ eV where — t. The dotted line shows how the spectrum would continue if Ts -C 10'' yr. The 
case of a bursting source corresponds to a slice of the image in the te — E plane, as indicated in Fig. ^ by 
dashed lines. For both spectra, d — 30Mpc, and 7 = 2.. 



time delay, rioo ^ 10 yr, as can be seen from the likelihood function marginalized over Ts 
and A^o (see Fig. [1^). According to Eq. (|T2|) this translates into a tentative bound for the 
r.m.s. magnetic field, namely. 



B <2x 10 
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-1/2 



d 



-1 



IMpc 



30Mpc 



(15) 



which also applies to magnetic fields in the halo of our Galaxy if d is replaced by the lesser 
of the source distance and the linear halo extent. If confirmed by future data, this bound 
would be at least two orders of magnitude more restrictive than the best existing bounds 
which come from Faraday rotation measurements [see Eq. (|14|)] and, for a homogeneous 



EGMF, from CMB anisotropics ||164|| . UHECRs are therefore at least as sensitive a probe 
of cosmic magnetic fields as other measures in the range near existing limits such as the 
polarization ||165|| and the small scale anisotropy ||166|] of the CMB. 



More generally, confirmation of a clustering of EHECRs would provide significant in- 
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Figure 10: The likelihood, £, marginalized over Tg and Nq as a function of the average time delay at 
10^° eV, Tioo, assuming a source distance d — 30Mpc. The panels are for pair # 3 through # 1, from top to 
bottom, of the AGASA pairs |148{ . Solid lines are for 7 = 1.5, dotted lines for 7 = 2.0, and dashed lines for 
7 = 2.5. 



formation on both the nature of the sources and on large-scale magnetic fields ||167|| . This 
has been shown quantitatively [|162|| by applying the hybrid Monte Carlo likelihood analysis 
discussed above to simulated clusters of a few tens of events as they would be expected 
from next generation experiments such as the High Resolution Fly's Eye the Tele- 
scope Array and most notably, the Pierre Auger Project (see Sect. 2), provided the 
clustering recently suggested by the AGASA experiment ||148| , |168|| is real. The proposed 
AirWatch type satellite observatory concepts ^ ^ might even allow one to detect 
clusters of hundreds of such events. 

classi- 



Five generic situations of UHECR time-energy images were discussed in Ref. 
fied according to the values of the time delay te induced by the magnetic field, the emission 
timescale of the source Ts, as compared to the lifetime of the experiment. The likelihood 
calculated for the simulated clusters in these cases presents different degeneracies between 
different parameters, which complicates the analysis. As an example, the likelihood is de- 
generate in the ratios Nq/Ts, or Nq/Atiqo, where Nq is the total fluence, and Arioo is the 
spread in arrival time; these ratios represent rates of detection. Another example is given 
by the degeneracy between the distance d and the injection energy spectrum index 7. Yet 
another is the ratio {(IteY^^ /h, that controls the size of the scatter around the mean of the 
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te — E correlation. Therefore, in most general cases, values for the different parameters 
cannot be pinned down, and generally, only domains of validity are found. In the following 
the reconstruction quality of the main parameters considered is summarized. 

The distance to the source can be obtained from the pion production signature, above 
the GZK cut-off, when the emission timescale of the source dominates over the time delay. 
Since the time delay decreases with increasing energy, the lower the energy Eq, defined by 
— Ts, the higher the accuracy on the distance d. The error on d is, in the best case, 
typically a factor 2, for one cluster of ~ 40 events. In this case, where the emission timescale 
dominates over the time delay at all observable energies, information on the magnetic field is 
only contained in the angular image, which was not systematically included in the likelihood 
analysis of Ref. ||162|| due to computational limits. Qualitatively, the size of the angular 



image is proportional to B{dlcy^^ / E, whereas the structure of the image, i.e., the number of 
separate images, is controlled by the ratio d^^'^B /l]./"^ / E. Finally, in the case when the time 
delay dominates over the emission timescale, with a time delay shorter than the lifetime of 
the experiment, one can also estimate the distance with reasonable accuracy. 

Some sensitivity to the injection spectrum index 7 exists whenever events are recorded 
over a sufficiently broad energy range. At least if the distance d is known, it is in general 
comparatively easy to rule out a hard injection spectrum if the actual 7 > 2.0, but much 
harder to distinguish between 7 = 2.0 and 2.5. 

If the lifetime of the experiment is the largest time scale involved, the strength of the 
magnetic field can only be obtained from the time-energy image because the angular image 
will not be resolvable. When the time delay dominates over the emission timescale, and is, 
at the same time, larger than the lifetime of the experiment, only a lower limit corresponding 
to this latter timescale, can be placed on the time delay and hence on the strength of the 
magnetic field. When combined with the Faraday rotation upper limit Eq. ([T^), this would 
nonetheless allow one to bracket the r.m.s. magnetic field strength within a few orders of 
magnitude. In this case also, significant information is contained in the angular image. If 
the emission time scale is larger then the delay time, the angular image is obviously the only 
source of information on the magnetic field strength. 

The coherence length l^. enters in the ratio {drEY^"^ /Ic that controls the scatter around the 
mean of the te — E correlation in the time-energy image. It can therefore be estimated from 
the width of this image, provided the emission timescale is much less than te (otherwise the 
correlation would not be seen), and some prior information on d and te is available. 

An emission timescale much larger than the experimental lifetime may be estimated if a 
lower cut-off in the spectrum is observable at an energy Eq, indicating that Ts ~ teq. The 
latter may, in turn, be estimated from the angular image size via Eq. ([T^), where the distance 
can be estimated from the spectrum visible above the GZK cut-off, as discussed above. An 



example of this scenario is shown in Fig. For angular resolutions A^, timescales in the 
range 

3x10^ {^X (^-\ yr <Ts~rs< 10^---10^ ( ^ \\i (16) 

V 1° ; vioMpc/-^ ~ ^ ^~ viooEevy ^ ^ ' 

could be probed. The lower limit follows from the requirement that it should be possible 
to estimate Te from Qe, using Eq. ([l^), otherwise only an upper limit on Tg, corresponding 
to this same number, would apply. The upper bound in Eq. (|16D comes from constraints 
on maximal time delays in cosmic magnetic fields, such as the Faraday rotation limit in 
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Figure 11: (a) Arrival time-energy histogram for 7 = 2.0, noo = 50 yr, Ts = 200 yr, — 1 Mpc, d = 50Mpc, 
corresponding to B ~ 3 x 10~^^ G. Contours are in steps of a factor 10" '' = 2.51; (b) Example of a cluster 
in the arrival time-energy plane resulting from the cut indicated in (a) by the dashed line at r ~ 100 yr; (c) 
The likelihood function, marginaHzed over Nq and 7, for d = 50 Mpc, Ic ~Mpc, for the cluster shown in (b), 
in the Tg — tiqo plane. The contours shown go from the maximum down to about 0.01 of the maximum in 
steps of a factor 10"'^ = 1.58. Note that the likelihood clearly favors Tg ~ 4tioo- For rioo large enough to 
be estimated from the angular image size, Tg ^ Tobs can, therefore, be estimated as well. 



the case of cosmological large-scale field (smaller number) and knowledge on stronger fields 
associated with the large-scale galaxy structure (larger number). Eq. (p!6D constitutes an 
interesting range of emission timescales for many conceivable scenarios of UHECRs. For 
example, the hot spots in certain powerful radio galaxies that have been suggested as UHECR 
sources ||169|| , have a size of only several kpc and could have an episodic activity on timescales 
of ~ 10^ yr. 

A detailed comparison of analytical estimates for the distributions of time delays, energies, 
and deflection angles of nucleons in weak random magnetic fields with the results of Monte 
Carlo simulations has been presented in Ref. [|1 70|| . In this work, deflection was simulated 



by solving a stochastic differential equation and observational consequences for the two 
major classes of source scenarios, namely continuous and impulsive UHECR production, 
were discussed. In agreement with earlier work ||150|| it was pointed out that at least in the 
impulsive production scenario and for an EGMF in the range 0.1 — 1 x 10^^ G, as required 
for cosmological GRB sources, there is a typical energy scale Eh ~ 10^°'^ — 10^^'^ eV below 
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which the flux is quasi-steady due to the spread in arrival times, whereas above which the 
flux is intermittent with only a few sources contributing. 



General Case 

Unfortunately, neither the diffusive limit nor the limit of nearly rectilinear propagation is 
likely to be applicable to the propagation of UHECRs around 10^° eV in general. This is 
because in magnetic fields in the range of a few 10~^ G, values that are realistic for the 
Supergalactic Plane ||157| , |156|| , the gyro radii of charged particles is of the order of a few 



Mpc which is comparable to the distance to the sources. An accurate, reliable treatment in 
this regime can only be achieved by numerical simulation. 

To this end, the Monte Carlo simulation approach of individual trajectories developed in 
Refs. [|161| , |162|] has recently been generalized to arbitrary deflections [|158|] . The Supergalactic 



Plane was modeled as a sheet with a thickness of a few Mpc and a Gaussian density profile. 
The same statistical description for the magnetic field was adopted as in Refs. ||161| , |162|| , 



but with a field power law index rin = —11/3, representing a turbulent Kolmogorov type 
spectrum, and weighted with the sheet density profile. It should be mentioned, however, 
that other spectra, such as the Kraichnan spectrum, corresponding to uh = —7/2, are also 
possible. The largest mode with non-zero power was taken to be the largest turbulent eddy 
whose size is roughly the sheet thickness. In addition, a coherent field component B^. is 
allowed that is parallel to the sheet and varies proportional to the density profile. 

When CR backreaction on the weakly turbulent magnetic field is neglected, the diffusion 
coefficient of CR of energy E is determined by the magnetic field power on wavelengths 
comparable to the particle Larmor radius, and can be approximated by 

As a consequence, for the Kolmogorov spectrum, in the diffusive regime, where te ^ d, 
the diffusion coefficient should scale with energy as D{E) oc E^/^ for Vg < L/{2ti), and as 
D(E) (X E in the so called Bohm diffusion regime, > L/(27r). This should be reflected 
in the dependence of the time delay te on energy E: From the rectilinear regime, te ^ d, 
hence at the largest energies, where te oc this should switch to te oc E~^ in the regime 
of Bohm diffusion, and eventually to te oc E~^^^ at the smallest energies, or largest time 



delays. Indeed, all three regimes can be seen in Fig. |I2] which shows an example of the 
distribution of arrival times and energies of UHECRs from a bursting source. 

The numerical results indicate an effective gyroradius that is roughly a factor 10 higher 
than the analytical estimate, with a correspondingly larger diffusion coefficient compared to 
Eq. (^). In addition, the fluctuations of the resulting spectra between different magnetic 
field realizations can be substantial. This is a result of the fact that most of the magnetic field 
power is on the largest scales where there are the fewest modes. These considerations mean 
that the applicability of analytical flux estimates of discrete sources in specific magnetic field 
configurations is rather limited. 

In a steady state situation, diffusion leads to a modification of the injection spectrum 
by roughly a factor te, at least in the absence of significant energy loss and for a homoge- 
neous, infinitely extended medium that can be described by a spatially constant diffusion 
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Figure 12: The distribution of time delays te and energies E for a burst with spectral index 7 = 2.4 at 
a distance d = lOMpc, similar to Fig. ||, but for the Supergalactic Plane scenario discussed in the text. 
The turbulent magnetic field component in the sheet center is -B = 3 x 10~^ G. Furthermore, a vanishing 
coherent field component is assumed. The inter-contour interval is 0.25 in the logarithm to base 10 of the 
distribution per logarithmic energy and time interval. The three regimes discussed in the text, te oc E~'^ in 



the rectilinear regime E > 200 EeV, te oc E ' 
Te oc E^^/^ for E < 60 EeV are clearly visible. 



in the Bohm diffusion regime 60 EeV < E < 200 EeV, and 



coefficient. Since in the non-diffusive regime the observed spectrum repeats the shape of the 
injection spectrum, a change to a flatter observed spectrum at high energies is expected in 
the transition region |171| . From the spectral point of view this suggests the possibility of 
explaining the observed UHECR flux above ~ 10 EeV including the highest energy events 



with only one discrete source [|172 



Angular images of discrete sources in a magnetized Supercluster in principle contain in- 
formation on the magnetic field structure. For the recently suggested field strengths between 
~ 10^® G and ~ lyuG the angular images are large enough to exploit that information with 
instruments of angular resolution in the degree range. An example where a transition from 
several images at low energies to one image at high energies allows one to estimate the 
magnetic field coherence scale is shown in Fig. |13. 
The newest AGASA data 
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however, indicate an isotropic distribution of EHECR. 
To explain this with only one discrete source would require the magnetic fields to be so 
strong that the flux beyond 10^° eV would most likely be too strongly suppressed by pion 
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production, as discussed above. The recent claim that the powerful radio galaxy Centaurus A 
at a distance of 3.4 Mpc from Earth can explain both observed flux and angular distributions 
of UHECRs above 10^'^'^eV for an extragalactic magnetic field of strength ~ 0.3/xG | 173 
was based on the diffusive approximation which is clearly shown by numerical simulations 
not to apply in this situation; the predicted distributions are too anisotropic even at 10^^ eV. 

This suggests a more continuous source distribution which may also still reproduce the 
observed UHECR flux above ~ 10^^ eV with only one spectral component |[L74|] . Statistical 



studies neglecting magnetic deflection can be performed analytically and also suggest that 
many sources should contribute to the UHECR spectrum ||175|| . A more systematic parameter 
study of sky maps and spectra in UHECR in different scenarios with magnetic fields is 
now being pursued ||176| , [L59|| . Such studies will also be needed to generalize analytical 
considerations on correlations of UHECR arrival directions with the large scale structure of 
galaxies | |1 77| | to scenarios involving significant magnetic deflection. 

Intriguingly, scenarios in which a diffuse source distribution follows the density in the 
Supergalactic Plane within a certain radius, can accomodate both the large scale isotropy 
(by diffusion) and the small scale clustering (by magnetic lensing) revealed by AGASA if a 
magnetic field of strength B > 0.5/iG permeates the Supercluster [|159 . 



Fig. in shows the distribution of arrival times and energies, the solid angle integrated 
spectrum, and the angular distribution of arrival directions in Galactic coordinates in such 
a scenario where the UHECR sources with spectral index 7 = 2.4 are distributed according 
to the matter density in the Local Supercluster, following a pancake profile with scale height 
of 5 Mpc and scale length 20 Mpc. The r.m.s. magnetic field has a Kolmogorov spectrum 
with a maximal field strength -Bmax = 5 x 10~^ G in the plane center, and also follows the 
matter density. The observer is within 2 Mpc of the Supergalactic Plane whose location 
is indicated by the solid line in the lower panel and at a distance d = 20 Mpc from the 
plane center. The absence of sources within 2 Mpc from the observer was assumed. The 
transition discussed above from the diffusive regime below ~ 2 x 10^° eV to the regime of 
almost rectilinear propagation above that energy is clearly visible. 

Detailed Monte Carlo simulations performed on these distributions reveal that the anisotropy 
decreases with increasing magnetic field strength due to diffusion and that small scale clus- 
tering increases with coherence and strength of the magnetic field due to magnetic lensing. 
Both anisotropy and clustering also increase with the (unknown) source distribution ra- 
dius. Furthermore, the discriminatory power between models with respect to anisotropy and 
clustering strongly increases with exposure ||159 . 



As a result, a diffuse source distribution associated with the Supergalactic Plane can 
explain most of the currently observed features of ultra-high energy cosmic rays at least 
for field strengths close to 0.5/iG. The large-scale anisotropy and the clustering predicted 
by this scenario will allow strong discrimination against other models with next generation 
experiments such as the Pierre Auger Project. 
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Figure 13: Angular image of a point-like source in a magnetized Supergalactic Plane, corresponding to one 
particular magnetic field realization with a maximal magnetic field in the plane center, Smax = 5 x 10~^ G, 
all other parameters being the same as in Fig. The image is shown in different energy ranges, as indicated, 
as seen by a detector of ~ f ° angular resolution. A transition from several images at lower energies to only 
one image at the highest energies occurs where the linear deflection becomes comparable to the effective field 
coherence length. The difference between neighboring shade levels is 0.1 in the logarithm to base 10 of the 
integral flux per solid angle. 
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Figure 14: The distribution of arrival times and energies (top), the soHd angle integrated spectrum (mid- 
dle, with 1 sigma error bars showing combined data from the Haverah Park [Q, the Fly's Eye and the 
AGASA experiments above 10^^ eV), and the angular distribution of arrival directions in Galactic coor- 
dinates (bottom, with color scale showing the intensity per solid angle) in the Supercluster scenario with 
continuous source distribution explained in the text, averaged over 4 magnetic field realizations with 20000 
particles each. 
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6 Conclusions 

Ultra-high energy cosmic rays have the potential to open a window to and act as probes of 
new particle physics beyond the Standard Model as well as processes occurring in the early 
Universe at energies close to the Grand Unification scale. Even if their origin will turn out 
to be attributable to astrophysical shock acceleration with no new physics involved, they 
will still be witnesses of one of the most energetic processes in the Universe. Furthermore, 
complementary to other methods such as Faraday rotation measurements, ultra-high energy 
cosmic rays can be used as probes of the poorly known large scale cosmic magnetic fields. 
The future appears promising and exciting due to the anticipated arrival of several large 
scale experiments. 
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